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Abstract. Josephson junctions have demonstrated enormous potential as qubits for scalable quan- 
tum computing architectures. Here we discuss the current approaches for making multi-qubit cir- 
cuits and performing quantum information processing with them. 
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P-i! 1. Introduction 
> 

^ . Macroscopic quantum behavior in a Josephson junction (JJ) was first demon- 
^ \ strated in the mid-1980's by John Clarke's group at UC Berkeley (Martinis et al., 
1985; Devoret et al., 1985; Martinis et al., 1987; Clarke et al., 1988). These exper- 
iments used a superconducting device referred to as a large area, current-biased 
JJ, which would later become the phase qubit. Beginning in the mid-1990's the 
group of James Lukens at SUNY Stony Brook (Rouse et al., 1995; Friedman 



. 5t , et al., 2000) and a collaboration between the Delft University group of Hans 

* mgeller@uga.edu 
epritchett@hal.physast.uga.edu 

* ats@math.uga.edu 

^ Present address: Physics Department and Insitute for Quantum Computing, University of 
Waterloo, Waterloo, Ontario N2L 3G1, Canada; fwilhelm@iqc . ca 




© 2008 Springer. Printed in the Netherlands. 



2 



GELLER ETAL. 



Mooij and the MIT group of Terry Orlando (Mooij et al., 1999; van der Wal 
at al., 2000) demonstrated macroscopic quantum behavior in superconducting 
loops interrupted by one or more JJs (called superconducting quantum interfer- 
ence devices, or SQUIDS), what would later become flux qubits. And in the 
late-1990's the group of Yasunobu Nakamura at NEC in Tsukuba (Nakamura 
et al., 1997; Nakamura et al., 1999) developed the first Cooper-pair box or charge 
qubit. Many of the earlier experiments were motivated by seminal theoretical work 
of Caldeira and Leggett (Caldeira and Leggett, 1981; Caldeira and Leggett, 1983). 

The modem era of superconducting quantum computation began in 2002. That 
year, the group of Siyuan Han at the University of Kansas and the group of John 
Martinis, then at NIST Boulder and currently at UC Santa Barbara, independently 
showed that long-lived quantum states in a current-biassed JJ can be control- 
lably prepared, manipulated, and subsequently measured (Yu et al., 2002; Martinis 
et al., 2002). This same year, the group of Michel Devoret, then at the CEA in 
Saclay and currently at Yale University, demonstrated similar quantum control 
using a Cooper-pair box (Vion et al., 2002). These experiments suggest that JJ- 
based qubits can be used as the building blocks of a solid-state quantum computer, 
creating a tremendous interest in this intrinsically scalable approach. An impres- 
sive list of additional experimental achievements soon followed, including the 
demonstration of two-qubit quantum logic (Yamamoto et al, 2003). 

In this chapter we will review the current approaches for making multi-qubit 
systems. For a more detailed discussion of single qubits we refer to the excellent 
review by MakhUn, Schon, and Shnirman (Makhlin et al., 2001). Also, a recent 
introductory account of the field has been given by You and Nori (You and Nori, 
2005). The approach we follow here is to construct circuit models for the basic 
qubits and coupled-qubit architectures. Many designs have been proposed, but 
only the simplest have been implemented experimentally to date. 

After reviewing in Sec. 2 the basic phase, flux, and charge qubits, we discuss 
three broad classes of coupling schemes. The simplest class uses fixed linear cou- 
pling elements, such as capacitors or inductors, and is discussed in Sec. 3. The 
principal eff'ect of fixed, weak couplings is to lift degeneracies of the uncoupled 
qubit pair. However, because such interactions are always present (always turned 
on), the uncoupled qubit states, which are often used as computational basis states, 
are not stationary. A variety of approaches have been proposed to overcome this 
shortcoming. In Sec. 4 we discuss tunable couplings that allow the interactions 
of Sec. 3 to be tuned, ideally between "on" and "off' values. A related class of 
dynamic coupUngs is discussed in Sec. 5, which make use of coupUng elements 
that themselves have active internal degrees of freedom. They act like tunable 
coupling elements, but also have additional functionahty coming from the abihty 
to excite the internal degrees of freedom. Examples of this are resonator-based 
couplings, which we discuss in some detail. 
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Figure 1. Circuit model for a current-biased JJ, neglecting dissipation. Here a = fi/2e. 



2. The basic qubits: phase, flux, and charge 



The primitive building block for all the qubits is the JJ shown in Fig. 1 . The low- 
energy dynamics of this system is governed by the phase difference (p between the 
condensate wave functions or order parameters on the two sides of the insulating 
barrier. The phase difference is an operator canonically conjugate to the Cooper- 
pair number difference A^, according to^ 



(1) 



The low-energy eigenstates if/m(<fi) of the JJ can be regarded as probability-amplitude 
distributions in <p. As will be explained below, the potential energy U {(p) of the JJ 
is manipulated by applying a bias current / to the junction, providing an external 
control of the quantum states i//m(f), including the qubit energy-level spacing Ae. 
The crossed box in Fig. 1 represents a "real" JJ. The cross alone represents a 
nonhnear element that satisfies the Josephson equations^ 



/ = /q sin ^ and V - a<p. 



(2) 



with critical current Iq. The capacitor accounts for junction charging.^ A single JJ 
is characterized by two energy scales, the Josephson coupling energy 



^ 2e 



(3) 



' We define the momentum P to be canonically conjugate to (p, and = P/ti. In the phase 

dip ' 



representation, N = -i4- 



2 a = hl2e. 

^ This provides a simple mean-field treatment of the wrer-condensate electron-electron interac- 
tion neglected in the standard tunneUng Hamiltonian formaUsm on which the Josephson equations 
are based. 
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where e is the magnitude of the electron charge, and the Cooper-pair charging 
energy 

£o - ^. (4) 

2C 

with C the junction capacitance. For example, 

320 neV 

£j = 2.05meVx7obuA] and E^ = -——, (5) 

C[pF] 

where /o[/uA] and C[pF] are the critical current and junction capacitance in mi- 
croamperes and picofarads, respectively. In the regimes of interest to quantum 
computation, and Ec are assumed to be larger than the thermal energy UqT but 
smaller than the superconducting energy gap Ajc, which is about ISOjueV in Al. 
The relative size of and E^. vary, depending on the specific qubit implementa- 
tion. 

2.1. PHASE QUBIT 

The basic phase qubit consists of a JJ with an external current bias, and is shown 
in Fig. 2. The classical Lagrangian for this circuit is 

1 o 

Ljj = -M^2 M = — . (6) 

Here ^ 

U - -Ej{ cos <p + s(p), with s = —, (7) 

^0 

is the efi'ective potential energy of the JJ, shown in Fig. 3. Note that the "mass" 
M in (6) actually has dimensions of mass x length'^. The form (6) results from 
equating the sum of the currents flowing through the capacitor and ideal Josephson 
element to /. The phase qubit implementation uses Ej » ^c. 
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Figure 3. Effective potential for a current-biased JJ. The slope of the cosine potential is s. The 
potential is harmonic for the qubit states unless s is very close to 1. 



According to tlie Joseplison equations, the classical canonical momentum P - 
|| is proportional to the charge Q or to the number of Cooper pairs Qjle on the 
capacitor according XoP - HQ/le. The quantum Hamiltonian can then be written 



as 



(8) 



where (p and A'^ are operators satisfying (1). Because U depends on s, which itself 
depends on time, Hjj is generally time-dependent. The low lying stationary states 
when s < I are shown in Fig. 4. The two lowest eigenstates |0) and |1) are used to 
make a qubit. Ae is the level spacing and AU is the height of the barrier. 

A useful "spin 5" form of the phase qubit Hamiltonian follows by projecting 
(8) to the qubit subspace. There are two natural ways of doing this. The first is to 
use the basis of the j-dependent eigenstates, in which case 



where 



(Or 



and a»po - sJlE^/H. 



(9) 



(10) 



The 5-dependent eigenstates are called instantaneous eigenstates, because s is usu- 
ally changing with time. The time-dependent Schrodinger equation in this basis 
contains additional terms coming from the time-dependence of the basis states 
themselves, which can be calculated in closed form in the harmonic limit (Geller 
and Cleland, 2005). These additional terms account for all nonadiabatic effects. 
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Figure 4. Effective potential in the anharmonic regime, with s very close to 1. State preparation 
and readout are carried out in this regime. 



The second spin form uses a basis of eigenstates with a fixed value of bias, sq. 
In this case 

H= '-^cr'-^{s-so)cr\ (11) 

where 

1 i2EcV 
^ = 4(1-^0)"^ and 4 = ■ (12) 

This form is restricted to |s - sqI 1> but it is very useful for describing rf pulses. 

The angle £ characterizes the width of the eigenstates in (p. For example, in the 
5o-eigenstate basis (and with in the harmonic regime), we have'* 

^ = ;toia""^ + arcsin(5o)cr°, with x„f,j' = {m\(p\m'). (13) 

Here x^m' is an effective dipole moment (with dimensions of angle, not length), 
and = {/ V2. 



2.2. CHARGE QUBIT 

In the charge qubit, the H current is provided capacitively, by changing the voltage 
Vg on a gate, as in Fig. 5. In this case Ej «; Ec, and the small capacitance is 
achieved by using a Cooper-pair box, which is a nanoscale superconducting island 
or quantum dot. 

The Lagrangian and Hamiltonian for this system are 

L ^ ^cn^iC + Cg)(p^ + Ej cos (f - aCgVg(p (14) 
cr" is the identity matrix. 
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Figure 5. Basic charge qubit circuit. The upper wire constitutes the superconducting box or island. 



and 



Here 



H = E^{N - N„f - Ej cos (f, with = ^^^^ 



2(C+Cg) 



2e 



(15) 



(16) 



is the gate charge, the charge qubit's control variable. 

It is most convenient to use the charge representation here, defined by the 
Cooper-pair number eigenstates |n) satisfying 



N\n) = n\n). 



(17) 



Because e"''|n) - \n + 1), the cos^ term in (15) acts as a Cooper-pair tunneling 
operator. In the qubit subspace, 



(A^ - Naf - (A^g - i)o^ + const. 



cos if 



(18) 
(19) 
(20) 



The charge qubit Hamiltonian can then be written in spin form in the {|0), |1)} 
charge basis as 



or in the |-)} basis of A^^g = ^ eigenstates 

10) ± |1) 



l±) 



(21) 



(22) 
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Figure 6. Basic rf-SQUID flux qubit and circuit model. <I> is the total flux threading the ring. The 
dashed curve in the upper figure indicates the integration contour T used to derive condition (26). 
The coil in the lower figure has self-inductance L. 



as 



(23) 



2.3. FLUX QUBIT 

The flux qubit uses states of quantized circulation, or magnetic flux, in a SQUID 
ring. The geometry is illustrated in Fig. 6. The current bias in this case is supplied 
by the circulating supercurrent. The total magnetic flux O can be written as 

0) = Ox - cLI, (24) 

where Ox is the external contribution and cLI is the self-induced component, with 

/ = aCip + Iq sin (p (25) 

the circulating current and L the self-inductance.^ The relations (24) and (25) 
determine O given (p, but there is a second condition relating these quantities, 

' L here is not to be confused with the Lagrangian. 
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Figure 7. Double-well potential of the flux qubit. The dashed curve is the cosine potential of the JJ 
alone; the solid curve shows the modification caused by the self-inductance of the ring. The states 
|0) and |1) are that of circulating and counter-circulating supercurrent, which become degenerate at 
the maximal frustration point (31). 



namely 



O w he 

— = mod 1, with = —. (26) 



This second condition follows from the Meissner effect, which says that the cur- 
rent density in the interior of the ring vanishes, requiring the total vector potential 
A to be proportional to the gradient of the phase of the local order parameter. It is 
obtained by integrating A around the contour T in Fig. 6. 
The relation (24) then becomes 

a'Cip + E,smcp+ - = 0, (27) 

where 

= (28) 

Vie 

This leads to the Lagrangian and Hamiltonian 

L.-a^C^^.E,cos,-^(,-^) (29) 

and 

H^E,N^-E,co.,.^[,-^). (30) 
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The ring's self-inductance has added a quadratic contribution to the potential 
energy, centered at 27rOx/Osc. 

The control variable in the flux qubit is Ox. By choosing 

Ox 1 

^.-modl, (31) 

one produces the double-well potential shown in Fig. 7. The condition (31) corre- 
sponds of the point of maximum frustration between the two directions of circu- 
lating supercurrent. By deviating slightly from the point (31), the energies of the 
|0) and |1) change, without changing the barrier height that controls the tunneling 
between the wells. 

We can write the flux qubit Hamiltonian in spin form as 

H = B^cr" + B^cr^, (32) 

where B^ and Bx are parameters that depend on the SQUID geometry and Ox. 
In the simplest rf SQUID flux qubit discussed here, B^ characterizes the well 
asymmetry, and is tunable (via Ox), whereas Bx depends on the barrier height 
and is fixed by the value of Ej. However, below we will describe a modification 
that allows the barrier height to be tuned as well. 

Hybrid charge-flux qubits have also been demonstrated, and have shown to be 
successful in reducing decoherence caused by interactions with the environment 
(Vion et al., 2002). 



3. Fixed linear couplings 



By fixed linear couplings we refer to coupling produced by electrically linear 
elements such as capacitors or inductors that lead to interaction Hamiltonians with 
fixed coupling strengths. In the cases usually considered, the coupling strengths 
are also weak, much smaller than the qubit level spacing, and we will assume that 
here as well. We discuss two prominent examples, capacitively coupled phase and 
charge qubits. For discussions of the third prominent example, inductively coupled 
flux qubits, we refer the reader to the literature (Mooij et al., 1999; Orlando et al., 
1999; MakhUn et al., 2001; Massen van den Brink, 2005). 



3.0.1. Capacitively coupled phase qubits 

Capacitively coupled phase qubits have been demonstrated by the University of 
Maryland group of Fred Wellstood (Berkley et al., 2003) and by the UC Santa 
Barbara group of John Martinis (McDermott et al., 2005). The architecture was 
discussed theoretically by Johnson et al. (Johnson et al, 2003), Blais et al. (Blais 
et al., 2003), and Strauch et al. (Strauch et al., 2003). 
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Figure 8. Capacitively coupled phase qubit circuit 



Referring to Fig. 8, the equations of motion for the two phase variables 

Qr^(Ci + Cint)ipi + SjiCsin ipi - si) - a^Cinm = 
a^(C2 + Cint)^2 + £'j2(sin (fii - si) - a^Qnt^i = 0, 

and the Lagrangian is 

^ - X [t^^; + Cint)(pf + Eji(cOS(pi + Si(pi)] - a^Cint(pi(p2- 



To find the Hamiltonian, invert the capacitance matrix in 

Pi \ _ 2 / Ci + C,nt -Cm \( <Pi 

Pi] *^ I -Cint C2 + Cint ]\(f2 



where the pi are the canonical momenta. This leads to 

.2 



i ' 



+ 



P\P2 



where 



Ci = C, + {Cr^,+C-^'r\ 

Qnt = CiC2(Cj ^ + + C^J) \ 



This can be written as 



^ a = n/2e 
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where 



{lef 



-int 



(42) 



The arrow in (42) applies to the further simpUfied case of identical qubits and 
weak coupling. 

The coupling constant g' defined in in (42) is inconvenient, however, because 
the energy scale appearing in (42) is too small. A better definition is 



g = 



Cint \^ 



where t is the scale introduced in (12). 

In the instantaneous basis, the spin form of the momentum operator is 



1 
-1 



where 



(43) 



(44) 



m = <0|p|l) = - 



Then 



(45) 



(46) 



In the uncoupled qubit basis {|00), |01), |10), |11)}, the qubit-qubit interaction in 
terms of (43) is simply 



SH^g 



^ 0-1 

10 

10 

-1 ; 



(47) 



Two-qubit quantum logic has not yet been demonstrated with this architec- 
ture. Methods for performing a controUed-Z and a modified swap gate have been 
proposed by Strauch et al. (Strauch et al., 2003), and four controUed-NOT imple- 
mentations have also been proposed recently (Geller et al., 2006). 



3.0.2. Capacitively coupled charge qubits 

A circuit for capacitively coupled charge qubits is given in Fig. 9. This architecture 
has been demonstrated by Pashkin etal. (Pashkin et al, 2003), and used to perform 
a CNOT by Yamamoto et al. (Yamamoto et al., 2003). This work is currently the 
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Figure 9. Capacitively coupled charge qubit circuit. 

most advanced in the field of solid-state quantum information processing. The 
equations of motion for the two phases are^ 

a^iCi + Cgi + Cint)<pi + En sin^i - aCgiVgi - a^Cimipi = (48) 
a^{C2 + Cgi + Cint)^2 + E}2 sin (p2 - aC^Vgi - a^Ci^m = 0, (49) 
and the Lagrangian is 

L^Yj [t(C'' + + Cint)^P^ + Eu cos - aCgiV^m] - ap-Cint(pm- (50) 



Then the Hamiltonian is 



\iPi + aCgiVgi)^ 



+ 



-E.oos^, 
(pi + aCg\Vgi){p2 + aCg2Vg2) 



a 



int 



where 



Cj = Ci + Cgi + 



C-'+iC2 + C„2) 



-1 



C2 = C2 + Cg2 + CrJ + (Ci + Cgi) 



Cint — 



Ci + Cgi + C2 + Cg2 + (Ci + Cgi)(C2 + Cg2)Cr„;. 
This can be written as 



where 



(J// = g(Afi - Afgi)(A^2 - %), £c,- = 



(f2)^ 
2C; 



g = 



(e2f 

2Cint 



(51) 

(52) 

(53) 
(54) 

(55) 
(56) 



a = n/2e 
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physical model circuit model 

Figure 10. Tuning £j with a dc SQUID. 
The spin form in the charge basis is 

i 

with 

Sti = f [(A^gi - 5)^2 + (^g2 - i)'^] + 1^X2- (58) 
When A^gi - Ng2 = ^, this is a pure Ising interaction. 

4. Tunable couplings 

By introducing more complicated coupUng elements, we can introduce some de- 
gree of tunability into the architectures discussed above. 

4.0.3. Tunable E] 

A simple way to make the Josephson energy E] effectively tunable in a circuit is 
to use a well known quantum interference effect in that occurs in a dc SQUID; see 
Fig. 10. The tunability of E] can be understood from two different viewpoints. 

The first is to imagine introducing a hole in a current-biased JJ as in the 
"physical" model of Fig. 10. Tunneling occurs in the up and down direction in 
each of the left and right arms of the interferometer. Recalling our interpretation 
of cos ^ as a Cooper-pair tunneling operator, the two arms of the interferometer 
result in 



COS (f 



+ e 



. *x 



-I- e 



-i{if- 



. *x 



(59) 
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Here we have assumed a symmetric interi'erometer. The first pair of terms corre- 
sponds to tunneling (in both the up and down directions) in the left arm, which 
acquires half of the total Aharonov-Bohm phase 2;rOx/Osc; the right arm has the 
opposite Aharonov-Bohm phase shift. Then the cos ip term in the potential energy 
of (8) becomes 



fij cos if £j(Ox) cos if, with £:j(Ox) = £j cos (^). (60) 

\ Csc ' 

The effective Josephson energy in (60) can be tuned by varying Ox. 

The second way to obtain (60) is to consider the circuit model in Fig. 10, and 
again assume symmetry (identical JJs). This leads to the coupled equations of 
motion 

aCip\ + /(jsin^i = 2 + /circ, (61) 
aClpi + /osin^2 = 5 -^irc- (62) 



Defining 



and using 



then leads to 



if = — - — (63) 



^i-^^^-oT ^^^^ 



a^(2C)ip + £'j(Ox) sin ip - al = 0, (65) 



in agreement with (60). 

The abihty to tune Ej is especially useful for inductively coupled flux qubits 
(Makhlin et al., 2001). 



4.0.4. Charge qubit register of Makhlin, Schon, and Shnirman 
Makhlin, Schon, and Shnirman have proposed coupling charge qubits by placing 
them in parallel with an inductor, such that the resulting LC oscillator (the capac- 
itance provided by the JJs) has a frequency much higher than the qubit frequency 
(Makhlin et al., 1999). The case of two qubits is illustrated in Fig. 11, but the 
method applies to more than two qubits as well. 

The derivation of the circuit Hamiltonian follows methods similar to that used 
above, and is 

-2 



LC ,E]iEj2 
+ ■* cr^V,. (66) 

40,2^2 12 
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Figure 12. Electrostatic transformer. 



The significant feature of the interaction in (66), compared to (58), is that the Ej's 
here can be tuned by using dc SQUIDs. This gives, in principle, a fully tunable 
interaction between any pair of qubits attached to the same inductor. 

4.0.5. Electrostatic transformer ofAverin andBruder 

Averin and Bruder (Averin and Bruder, 2003) considered a related coupled charge 
qubit circuit, shown in Fig. 12, which we have reorganized to emphasize the 
similarity to Fig. 1 1 . The Hamiltonian in this case is 



+ ^trans (67) 



and 



^trans = E^(N -N„- qf - £j COS (p, (68) 
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where 

q = q^-iNi- ^gi +N2- (69) 

% ^ 2A^g(i - (70) 

The operator q here is a function of the charge qubit variables, but commutes with 
the transformer degrees of freedom. 

As in the register of Makhhn, Schon, and Shnirman, we assume the trans- 
former degrees of freedom are fast compared with the qubit variables, so that the 
transformer remains in its instantaneous ground state manifold. Then 

-f^trans ^ ^oiq). (71) 

This finally leads to an effective Hamiltonian 



(72) 



involving charge qubit variables only, where 



a ^ (73) 

b = ^ ^ . (74) 

The discrete second-order derivative b, which can be interpreted as a capacitance, 
can be tuned to zero by varying qQ, providing the desired tunabiUty. 



4.1. RF COUPLING 



Finally, we briefly mention an interesting proposal by Rigetti, Blais, and Devoret, 
to use rf pulses to effectively bring permanently detuned qubits into resonance 
(Rigetti et al., 2005). This is a very promising approach, but has not yet been 
demonstrated experimentally. 



5. Dynamic couplings: Resonator coupled qubits 



Several investigators have proposed the use of LC resonators (Shnirman et al., 
1997; Makhlin et al., 1999; Mooij et al, 1999; You et al, 2002; Yukon, 2002; Blais 
et al., 2003; Plastina and Falci, 2003; Zhou et al., 2004), superconducting cavities 
(Blais et al, 2004; Wallraff et al, 2004), or other types of oscillators (Marquardt 
and Bruder, 2001; Zhu et al., 2003) to couple JJs together. Although harmonic 
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Figure 13. Two current-biased Josephson junctions (crossed boxes) coupled to a piezoelectric disc 
resonator. 

oscillators are ineffective as computational qubits, because the lowest pair of 
levels cannot be frequency selected by an external driving field, they are quite 
desirable as bus qubits or coupling elements. Resonators provide for additional 
functionality in the coupling, and can be made to have very high Q factor. Here 
we will focus on phase qubits coupled by nanomechanical resonators (Cleland 
and Geller, 2004; Geller and Cleland, 2005; Somborger et al., 2004; Pritchett and 
Geller, 2005). 

5.1. QUBIT-RESONATOR HAMILTONIAN 

The Hamiltonian that describes the low-energy dynamics of a single large-area, 
current-biased JJ, coupled to a piezoelectric nanoelectromechanical disk resonator, 
can be written as (Cleland and Geller, 2004; Geller and Cleland, 2005) 

^ = ^"'^mCm + fiajQa^a - ig ^ x^rn'c'lcm'ia - a'*'), (75) 

m mm' 

where the {cl,} and {c„,) denote particle creation and annihilation operators for 
the Josephson junction states (m = 0, 1,2, . . .),a and denote ladder operators 
for the phonon states of the resonator's dilatational (thickness oscillation) mode 
of frequency ajQ, g is a coupling constant with dimensions of energy, and x^m' = 
{m\ip\m'). The value of g depends on material properties and size of the resonator, 
and can be designed to achieve a wide range of values. An illustration showing 
two phase qubits coupled to the same resonator is given in Fig. 13. Interactions 
between the JJ and resonator may be controlled by changing the JJ current, giving 
rise to changes in the JJ energy spacing. At. For instance, a state can be transferred 
from the JJ to the resonator by bringing the JJ and resonator in resonance, Ae = 
fiojQ, and waiting for a specified period. 
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5.2. STRONG COUPLING AND THE RWA 

For small couplings g <«c Ae, the JJ-resonator system may be approximated by the 
Jaynes-Cummings model; this is usually referred to as the rotating wave approx- 
imation (RWA). However, once the coupling becomes comparable to the level 
spacing, g x Ae, the RWA breaks down. When the JJ is weakly coupled to the 
resonator, with g/Ae below a few percent, gates such as a memory operation 
(state transfer to and from the resonator) work well, and qubits are stored and 
retrieved with high fidelity. However, such gates are intrinsically slow. As g/Ae is 
increased, making the gate faster, the fidelity becomes very poor, and it becomes 
necessary to deviate from the RWA protocol. Below, we first discuss an analytical 
approach to capture the leading corrections to the RWA at intermediate coupling 
strengths (Sornborger et al., 2004). We then discuss a strong coupling information 
processing example: a quantum memory register (Pritchett and Geller, 2005). 

5.3. BEYOND THE RWA 

For simplicity we will consider only two levels in a single junction. However, 

all possible phonon-number states are included. The Hamiltonian may then be 
written as the sum of two terms, H - Hjc + V. The first term, 

Hjc = €0 cjco + ei cjci + Hojo a' a - igxm {c\cQa - CpCia^], (76) 

is the exactly solvable Jaynes-Cummings Hamiltonian, the eigenfunctions of which 
are known as dressed states. We will consider the second term, 

V = -jg[xoocJco(a - a") + xoicjcia - XQ\c\cQa' + xnc\c\{a - a')], (77) 

as a perturbation. The RWA applied to the Hamiltonian H amounts to neglect- 
ing V. Therefore, perturbatively including V is equivalent to perturbatively going 
beyond the RWA. 

5.3.1. Dressed states 

The eigenstates of Hjc, or the dressed states, are labeled by the nonnegative 
integers j - 0,1,2, . . . and a sign cr = ±1. On resonance, these are 

i,(rx _ io,j + i)-io-|i.y) , „, 

\>Pj} = ^ , = 0) (78) 

and 

Wj = eo + ij+ mojo + a yfJTT (a>d = 0) (79) 

Here, the vacuum (j-0) Rabi frequency on resonance is Qo(0) - 2g|xoi|/^. 
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5.3.2. Dressed state propagator 

In quantum computing applications one will often be interested in calculating 
transition amplitudes of the form 

{i\e-'""%), (80) 

where |i) and |f ) are arbitrary initial and final states of the uncoupled qubit-resonator 
system. Expanding |i) and |f ) in the dressed-state basis reduces the time-evolution 
problem to that of calculating the quantity 

C'.j'it) ^ {ip^le-'"'"'],^^:}, (81) 

as well as <^J|e-'^'/^|00) and (00|e-''^'/^|00). G'jfit) is a propagator in the dressed- 
state basis, and would be equal to da-a-'Sjfe''^' '^^ if V were absent, that is, in the 
RWA. 

To be specific, we imagine preparing the system at f = in the state |10), 
which corresponds to the qubit in the excited state m = 1 and the resonator in the 
ground state n - 0. We then calculate the interaction-representation probabihty 
ampUtude 

Cnnit) = e'^""^"^(mn\e-'""^m (82) 

for the system at a later time t to be in the state \mn). Here = + nhojQ. 
Inserting complete sets of the dressed states leads to 

emit) = ^<^/.J|I0)<00k-'^^/^|^/.p, (83) 
0-7 

and, for mn 00, 

2. \ (,r,\mn) ] \ G- G-- ] \ <^-|10) I " ^^^^ 



So far everything is exact within the model defined in Eq. (75). 

To proceed, we expand the dressed-state propagator in a basis of exact eigen- 
states \^a) of H, leading to 

G'^jfit) = Y^rj\'Va) l^a)* e-^'-'l'^. (85) 

a 

Here is the energy of stationary state \^a)- The propagator is an infinite sum of 
periodic functions of time. We approximate this quantity by evaluating the |*Pq.) 
and £q, perturbatively in the dressed-state basis. 

We test our perturbed dressed-state method for the case of a finite-dimensional 
single-qubit, five-phonon system. The bias current is chosen to make the system 
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Figure 14. Probabilities |cioP and |coiP for the intermediate case of ^/Ae = 0.30. Here there are 
large deviations from the RWA behavior, which are correctly accounted for by the dressed-state 
perturbative method. 

exactly in resonance. The Hamiltonian for this system is diagonaUzed numeri- 
cally, and the probability amplitudes c,„,i{t) are calculated exactly, providing a test 
of the accuracy of the analytic perturbative solutions. Setting the initial state to be 
Cmn(O) = Smi6nQ, as assumed previously, we simulate the transfer of a qubit from 
the Josephson junction to the resonator, by leaving the systems in resonance for 
half a vacuum Rabi period nfilg\xQ\\. 

In Fig. 14, we plot the probabilities for a relatively strong coupling, g/Ae = 
0.30. For this coupUng strength, the RWA is observed to fail. For example, the 
RWA predicts a perfect state transfer between the junction and the resonator, and 
does not exhibit the oscillations present in the exact solution. The dressed state 
perturbation theory does correctly capture these oscillations. 

5 .4. MEMORY OPERATION WITH STRONG COUPLING 

Here we study a complete memory operation, where the qubit is stored in the 
resonator and then transferred back to the JJ, for a large range of JJ-resonator 
couphng strengths (Pritchett and Geller, 2005). Also, we show that a dramatic 
improvement in memory performance can be obtained by a numerical optimiza- 
tion procedure where the resonant interaction times and off-resonant detunings are 
varied to maximize the overall gate fidelity. This allows larger JJ-resonator cou- 
plings to be used, leading to faster gates and therefore more operations carried out 
within the available coherence time. The results suggest that it should be possible 
to demonstrate a fast quantum memory using existing superconducting circuits, 
which would be a significant accomphshment in sohd-state quantum computation. 
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dimensionless interaction strength (in %) 

Figure 15. (upper panel) Memory fidelity for equator state 2"J(|0> + |1)) as a function of g/Ae, 
using botii tlie RWA (unfilled circles) and optimized (solid circles) pulse times, (lower panel) Time 
needed to store and retrieve state, using both the RWA (dashed curve) and optimized (solid curve) 
pulse times. 



In the upper panel of Fig. 15 we plot the memory fideUty for the qubit state 
2~2(|0) + |1)) as a function of g/Ae. We actually report the fidelity squared, 

F2 = |a*coo(?f)+y8*cio(?f)|' (86) 

which is the probability that the memory device operates correctly. As expected, 
the fidelity gradually decreases with increasing g. The lower panel of Fig. 15 gives 
the gate time as a function of g/Ae. These results suggest that memory fidelities 
better than 90% can be achieved using phase qubits and resonators with coherence 
times longer than a few tens of ns. 
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